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Abstract. We prove that the spectrum of an n-dimensional semiclassical radial Schrodinger 
operator determines the potential within a large class of potentials for which we assume no 
symmetry or analyticity. Our proof is based on the first two semiclassical trace invariants 
and on the isoperimetric inequality. 



1. Introduction 

In this paper we study the problem of determining a potential V in the class 

{V G C°°(R n ) : -h 2 A + V is selfadjoint on a domain containing C °°(M n )} (1.1) 

from its eigenvalues. Note that this class includes the set of potentials V which are bounded 
from below, because in this case the Friedrichs extension satisfies the condition. 

We prove that radial, monotonic potentials are spectrally determined within this class. 

Let n > 2. Let V G C°°(lR n ) be radial near in the sense that V (x) = R(\x\) for |x| < R 
for some Ro > 0, and suppose R satisfies R(0) = and R'{r) > for r G (0, Rq). Let 
Ao = R(Rq), and suppose Vo > Ao for \x\ > Rq. 

Theorem. IfV is in the class ( II. ip and the spectrum of — h 2 A + V agrees with the spectrum 
of —h 2 A + Vo in the interval (—00, Ao) up to order o{h 2 ) for h G (0, ho] for some h > 0, 
then there exists xo G M™ such that V(x — Xq) = Vo(x) for \x\ < Ro- 

In particular, the harmonic oscillator is spectrally determined. 

Corollary. If V is in the class (11. ip and if the spectrum of —h 2 A + V agrees with the 
spectrum of —h 2 A + \x\ 2 up to order o(h 2 ) for h G (0, ho] for some ho > 0, then there exists 
Xq G W 1 such that V(x — x ) = \x\ 2 . 



Date: February 11, 2011. 

Key words and phrases, inverse spectral theory, trace invariants, semiclassical Schrodinger operators. 

The first author is partially supported by a National Science Foundation postdoctoral fellowship, and the 
second author is partially supported by the National Science Foundation under grant DMS-0969745. The 
authors are grateful for the hospitality of the Mathematical Sciences Research Institute, where part of this 
research was carried out. 



2 



KIRIL DATCHEV, HAMID HEZARI, AND IVAN VENTURA 



The oldest inverse spectral result is that balls in R n are determined by their Dirichlet or 
Neumann eigenvalues within the class of bounded domains with smooth boundary. This 
follows from the isoperimetric inequality and the fact that the first two heat invariants, i.e. 
the first two terms of the asymptotic expansion of Tre iA as t ->■ 0+ give the volume of 
the domain and the area of its boundary. All other positive results for spectral uniqueness 
of domains rely on analyticity or symmetry assumptions on the domain and on the class of 
domains under consideration. In |Zel09j , Zelditch shows that any analytic domain in M 2 with 
one reflection symmetry is determined within the class of all such domains, and in |HeZelO] , 
the second author and Zelditch generalize this result to higher dimensions. In |HeZel 0p|, 
the second author and Zelditch show that an isospectral deformation of an ellipse within the 
class of smooth domains with two reflection symmetries is infinitesimally trivial. 

In general it is impossible to determine a potential V from the spectrum of the nonsemi- 
classical Schrodinger operator —A + V. For example, in |McTr81j . McKean and Trubowitz 
find an infinite dimensional family of potentials in C°°(R) which are isospectral with the 
harmonic oscillator V(x) = x 2 . However, if we have information about the spectrum of the 
semiclassical Schrodinger operator —h 2 A + V for all h in (0, ho], we can say more about the 
potential. In fact, in the Theorem above, it is sufficient for example to know the spectrum up 
to o{h 2 ) for h e {hj}™L 1 for a sequence hj — > 0. In a physical interpretation of the problem, 
varying h corresponds to varying a parameter of the system, such as a mass or coupling 
constant. 

The semiclassical Schrodinger operator is, in the following sense, the natural spectral ana- 
logue of the Laplacian on a domain or manifold. A semiclassical parameter h can be intro- 
duced into any spectral problem about the Laplacian on a domain or manifold by scaling 
the metric by h~ 2 , and hence producing a factor of h 2 in front of the Laplacian. Because of 
the homogeneity of the symbol of the Laplacian in the fibers of the cotangent bundle, the 
spectrum of the semiclassical Laplacian for any value of h can be deduced from the spectrum 
at h — 1: 

— Auj = XjUj -<=>- —h 2 Auj = Ej(h)uj, for all h > 0, where Ej(h) = h 2 \j. (1-2) 

The symbol of a Schrodinger operator is not homogeneous in the fibers, and so in this setting 
there is no analogue of (11. 2p . This partially accounts for the large family of isospectral 
potentials in the example of McKean and Trubowitz mentioned above: the nonsemiclassical 
spectrum alone provides much less information in this setting. 

In |IaSj Zw02] , Iantchenko, Sjostrand and Zworski study spectral inverse problems of semi- 
classical Schrodinger operators, recovering in particular Birkhoff normal forms of closed 
orbits. In |GuUr 07j. Guillemin and Uribe show that any real analytic potential V, symmet- 
ric with respect to all coordinate axes and with a unique global minimum, is determined 
within the class of all such potentials by the spectrum at the bottom of the well of its semi- 
classical Schrodinger operator. In [Hcz09], the second author removes one of the symmetry 
assumptions, giving a result with no symmetry assumptions in dimension 1, but keeps the 
assumption of analyticity. In |CoGu08p|, Colin de Verdiere and Guillemin give a differ- 
ent proof of the one dimensional result. In |Col0 8p|, Colin de Verdiere removes both the 
symmetry and analyticity assumptions in the one dimensional case, but adds a genericity 
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assumption. In |GuWa09p1 (see also |GuStllp"| §10.6]) Guillemin and Wang prove a one 



dimensional version of the result in the present paper: using the two trace invariants 

/ dxd£, [ \VV(x)\ 2 dxd£, (1.3) 

they show that an even function (or a suitable noneven function) is determined by its spec- 
trum within the class of functions monotonic away from 0. To the authors' knowledge, the 
result in the present paper is the first positive result for a semiclassical Schrodinger operator 
in higher dimensions with no analyticity, symmetry or genericity assumption on the class 
within which the potential is spectrally determined. 

We mention also some related results for determining an obstacle from its resonances. In 
[Ha Zw99j . Hassell and Zworski show that the 3-ball as an obstacle is determined by its 
Dirichlet resonances. In [C hr08 ] , Christiansen extends the result to n-balls and Neumann 
resonances. 



Our proof is based on the isoperimetric inequality and on the first two terms of a semiclassical 
trace formula due to Helffer and Robert |HeRo83] . In this it is very similar to the proof of 
spectral uniqueness of balls in M™, which uses the isoperimetric inequality and the first two 
terms of the heat trace. Because we use only the first two terms of the expansion, we only 
need to know the spectrum up to o{h 2 ). 

We are grateful to Maciej Zworski for his encouragement, and for several interesting dis- 
cussions about this paper and about the motivation for the inverse spectral problem for 
semiclassical Schrodinger operators. 



2. Review of semiclassical trace formulas 



For V G C°°(R n ; [0,oo)), let 

P d £ _/,2 A + y 

If P has discrete spectrum in [0,Ao), then for each h > 0, f(P) (defined by the functional 
calculus for self adjoint operators) has finite rank for every / 6 Co°((— oo, Ao)), and the trace 
of f(P) is determined by the eigenvalues of P. In this section we review three methods of 
proof for the trace formula 

Tr(/(P)) = 

7^(7 fm 2 + v)dxd^+^- ! \vv\ 2 f^\\^+v)dxd^+o(h^). {2A) 

{27Th) n \J R 2n 12 J R 2n J 

The two invariants in (jl.3p are obtained from this formula by approximating the character- 
istic function of an interval in C£°(R). 



The trace formula (12. ip comes from an asymptotic expansion of a(f(P)), the full semiclassical 
symbol of f(P) (here the Weyl symbol, or the standard symbol, or another symbol may be 
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used), of the form 

a(f(P)) = a (f(P)) + hcrxUiP)) + h 2 a 2 (f(P)) + ■■■ (2.2) 

combined with a version of Lidskii's formula for the trace of an operator as the integral along 
the diagonal of its integral kernel, as in for example |DiSj99, Theorem 9.6]: 



Tr(/(P)) = — ^ (7 M/(P)) + ha 1 {f(P)) + h 2 a 2 (f(P)) + ■ ■ ■} dxd^j . 

The terms of the expansion (12. 2p were first computed by Helffer and Robert |HeRo83l Propo- 
sition 5.3] using the Mellin inversion formula 



I(P) = ^ I /•'(.«) P-V.s. 

' (7 — 200 



2ni 



where F is the Mellin transform of /, given by F(s) = f£°t s l f{t)dt. In |DiSj99 The- 



orem 8.7] Dimassi and Sjostrand obtain an expansion using the Helffer- Sjostr and formula 
( |HeSj89| , see also for example |DiSj99[ (8.3)]) 



f(P) = -- [ df{z){z-P)- l dxd yi 
K Jc 



where / is an almost analytic extension of /. More recently, in |GuStllp §10], Guillemin 
and Sternberg use the Fourier inversion formula 



f(P) = i- / f(t)e itp dt 

27T J R 



to obtain the same expansion. This last method is a semiclassical version of Taylor's construc- 
tion in |Tay81 §12.1]. In each case this reduces the problem to computing the asymptotics 



of the symbol of one of P s , (z — P) 1 or e ltp (in fact in the first case Helffer and Robert 
use the contour integral formula P~ s = (27ri) _1 J z~ s (z — P)~ l dz for suitable 7 to reduce 
the problem to studying (z — P)^ 1 )- Using one of these methods, and either the standard or 
Weyl symbol, we obtain an expansion of the form (12. 2p . Regardless of the choice of symbol, 
we have a (f(P)) = f(a (P)), i.e. the principal symbol of f(P) is / of the principal symbol 
of P. Using the formulas for a ,ai,a 2 ,cr 3 derived in |GuStllp] §10.5] gives (12.11) . 



3. Proof of Theorem 



The spectrum of —h 2 A + Vq is contained in (0, 00) and is discrete in [0,Ao). Since the 
spectrum of —h 2 A + V is the same as the spectrum of —h 2 A + V , up to order o{h 2 ) on 
(— 00, A ), we conclude that < — h 2 A + V + o{h 2 ). It follows that V > 0. Indeed, otherwise 
we could take a function u G C^°(M ri ) with suppw C {V < 0}, giving 

0< f(-h 2 A + V)\u\ 2 dx + o(h 2 )\\u\\ 2 L2 = I V>| 2 + 0(/i 2 )|H|^i < 0, 

J J supp u 

for h sufficiently small, which is a contradiction. 
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Because the spectrum of —h 2 A + V is known on (—00, Ao) up to order o(h 2 ) we obtain from 
( 12. ip the two trace invariants in (11. 3p for each A G (0, Ao). By integrating in the £ variable, 
we rewrite these invariants as follows: 

(X-V) n/2 dx } [ \VV(x)\ 2 (X-V) n/2 dx. (3.1) 

{V{x)<\} J{V(x)<\} 

From the first invariant of (13. ip it follows that 

{x: V(x) < X,W(x) = 0} has Lebesgue measure zero. (3.2) 
We will prove this in §H 



Using (13.21) and the coarea formula we rewrite the invariants in (13.11) as 

r x ( r (a - vy' 2 \ ^ 

JO \J{V=s,X7V^0} 



iwi 



-dS^jds, (^J |VU|(A - V) n/2 dS^j ds. 



Using the fact that V = s in the inner integrand, the factor of (A — V) n / 2 = (A — s) n / 2 can 
be taken out of the surface integral, leaving 

(A- s) n / 2 h{s)ds, / (A - s) n/2 I 2 (s)ds, (3.3) 



Jo 



where 

h(s) = [ T^77\dS, I 2 (s) = [ \VV\dS. (3.4) 

J{V=s,VV^0} I Vl/ I J{V=s} 

We denote the integrals ( 13.31) by Ai +ri / 2 (/i)(A) and Ai +n / 2 (Ji)(A). These are Abel fractional 
integrals of 1\ and I2 (see for example [Zel98[ §5.2] and |GuStllp"l , (10.45)]), and they can 
be inverted by applying Ai+ n /2, using the formula 

-A a o r-^rAp = — 1 a . A a+ p, (3.5) 



r(a) a v(/3) p r(a + py 

and differentiating n + 1 times. From this we conclude that the functions I\ and I2 in ( 13.41) 
are spectral invariants for every s G (0, Ao). 

Integrating Ii and using the coarea formula again as well as ( 13.21) we find that the volumes 
of the sets {V < s} are spectral invariants for each s G (0, A ): 

I^ds'^ [ [ -r^—dSds' = [ ldx. (3.6) 

Jo J{V=s'yV^0} |VK I J{V<s} 

From Cauchy-Schwarz and the fact that 1% and I2 are spectral invariants we obtain 

([ Ids] < [ -^—dsf \VV\dS= [ -LdS [ R'dS, (3.7) 

\J{V=s} J J{V=s} I V 7 1/ I J{v=s} J{R=s} R J{R=s} 



6 KIRIL DATCHEV, HAMID HEZARI, AND IVAN VENTURA 

when s is not a critical value of V, i.e. by Sard's theorem for almost every s G (0, Ao). On 
the other hand, using the invariants obtained in (13. 6p and the fact that the sets {R < s} are 
balls, by the isoperimetric inequality (see for example |Oss78t pi 188]) we find 

ldS< [ IdS. 

{R=s} J{V=s} 

However, 

IdS) = [ -LdS ( R'dS, 



{R=s} J J{R=s} R J{R=s} 

because 1/R' and R' are constant on {R = s}. Consequently 

ldS= [ IdS, 

'{R=s} J{V=s} 

and so {V = s} is a sphere for almost every s, because only spheres extremize the isoperi- 
metric inequality. Moreover, 

2 

IdS I / -Jr-dS I \YV\dS. 



{V=s} / J{V=s} |VV| J{V=s} 

and so |VK| _1 and |W| are proportional on the surface {V = s} for almost every s, again 
by Cauchy-Schwarz. Using the equation (13. 7p to determine the constant of proportionality, 
we find that 

|W| 2 = R'iR-^s)) 2 = (R-^'is)- 2 = F(s) 
on {V = s}. In other words 

|VVf = F(V), (3.8) 

for all x G y _1 (s) for almost all s G (0,A ). However, because F(V) ^ when V ^ 0, it 
follows by continuity that this equation holds for all x G U _1 ((0, Ao)). 

To solve this equation, we restrict it to flowlines of W, with initial conditions taken on a 
fixed level set {V = s }: 

'x(t) = VV(x(t)), 
x(0) G {V = s }. 

Observe first that this flow takes level sets to level sets. Indeed, 

jV{x{t)) = W(x(t)) ■ x(t) = \VV(x(t))\ 2 = F(V(x(t))), 



(3.9) 



and this equation can be solved: 



v dV> 



i'M= I jjyry v(x(t)) = r\t). 



Note that / depends on s but not on x(0). 

We now use this to show that the integral curves x(t) are lines. Differentiating (I3.9P we find 

x {t) = ^sjV{x{t)) = (Hess ■ VV(x(t)). 
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But differentiating (13. 8 p gives 

2(HessV)- VV = F'(V)VV, 



so that 




Integrating this equation we find the integral curves are lines: 




VF(I^(t)) _ ( ) VF(V(x(t))) 
\x(0)\ H ' \x(0)\ 



Moreover, these lines are all directed at the center of the sphere {V = s Q } and have a speed 
dependent only on s Q but not on x(0), and consequently all level sets in the flowout of 
{V = s } are spheres with the same center. Since VV ^ on {0 < V < A } thanks to (13.8 j) . 
the integral curves cover this region, and we conclude that, up to a translation, V is radial 



on {\x\ < Ro}. Since the volumes (I3.6P are spectral invariants, it follows that V(x) = R(\x\) 
on {\x\ < R }. 



In this section we prove (I3.2p . To do this we will use the first invariant of (13.11) to show that 
V*dx, the pushforward of Lebesgue measure by V, is absolutely continuous with respect to 
Lebesgue measure on (— oo, Ao). From this it follows that V~ 1 {E) is Lebesgue-null whenever 
E C (— oo, Ao) is Lebesgue-null. Since the set of critical values of V is Lebesgue-null thanks 
to Sard's theorem, this completes the proof. 

To show that V*dx is absolutely continuous, we write the first invariant of (13. ip as 



for every A G (0, Ao), where xe denotes the characteristic function of E. The finiteness of the 
integrals in (14. ip implies that V*dx(— oo, A) < oo for each A < Ao, and hence V*dx is a finite, 
regular, Borel measure on each (— oo, A). Applying the Abel fractional integral transform as 
in fj3j we find that 



is a spectral invariant for every A G (0, Ao). For this we used (13.51) . which is proved using only 
the Fubini-Tonelli theorem followed by an elementary single-variable integral computation, 
and hence is valid here. 

We now show that, for m G {1, . . . , n}, 



4. Analysis of the critical set 




(4.1) 




d 

dX 



/(A - s) m X[ o,x)(s)(V*dx) = m /(A - s^xioMiV.dx), 




(4.2) 
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where the derivative is guaranteed to exist since, by induction, these are each spectral invari- 
ants. To prove ( 14. 2p . we write 



(X + h— s) m X[ a,x+h)(s)(V*dx) - / (A — s) m X[ a,x)(s)(V*dx) = 

JR 

[(X + h — s) m - (A - s) m ] X [o,x +h )(s)(V*dx) + /(A - s) m [ X[ o,x +h) (s) - X [o,x)(s)] (V.dx). 



For the first term we have 

\m— 1„ 



[(X + h- s) m - (A - s) m ] X [ ,x + h)( s )( V * dx ) =mh (X- s) m - 1 x M+fc )(s)(Ktfe) + O(^), 

thanks to the fact that J X[o,x+h)(s) (V*dx) < oo for small h. Using the monotone convergence 
theorem we find that 

[J(X-s) m x X[o,x)W(V m dx), h^O . 

Meanwhile for the second term we have, 

m w r / \ nl/I ,,> (/(A - s) m X{A}(s)(K^) = 0, h -»■ + , 

(A-s) m [x[o,A+^)( s )-X[o,A)(^)] (Kdx)^- ^ v 1 n ^^ _ 

Since the derivative in ( 14. 2 j) must exist, the limits from the left and right in ( 14. 3 p must be 
equal, and ( 14. 2 p is proved. In particular, we find that 



X[o,x)(s)(V*dx) = / X[o,A)(s)((V )*(ix) 
Jr 

for every A > 0, implying that V*dx agrees with an absolutely continuous measure on a 
family of sets which is generating for the Borel subsets of (— oo, Ao), and in particular that 
the critical set of V is Lebesgue-null. 
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